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Summary  of  2012-2013  Progress 
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•  self-consistent  simulations  of  pulsed  nanosecond  discharges  in  air. 

•  detailed  validation  with  experiments  and  analytical  model  results. 

•demonstration  of  volumetric  plasma  heating  and  radical  production  of  critical 
importance  in  combustion  applications. 

Ignition  of  H2-,  CH4-,  and  C2H4-Air  Mixtures 
•critical  assessment  of  plasma  kinetic  models  through  comparison  of  OH  decay  rates 
after  a  burst  of  nanosecond  pulses  below  ignition  threshold  temperatures  (~500  K). 
•importance  of  local  plasma  chemistry  effects  over  heat  transport  in  achieving 
“volumetric”  ignition  using  pulse  nanosecond  discharges. 

•detailed  parametric  studies  on  the  sensitivity  of  nanosecond  plasma  ignition  to  pressure, 
eq.  ratio,  pulsing  frequency,  burst  size,  initial  temperature,  and  dielectric  properties. 

Ignition  of  Heavy  Fuels  (n-Heptane) 

•  effect  of  nanosecond  plasma  on  the  two-stage  n-heptane  ignition  process. 
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Plasma-Coupled  Premixed  Flames 

construction  of  plasma  flame  kinetic  mechanisms,  including  electron  impact  dynamics 
of  all  major  species  in  flame  environments  (both  reactants  and  products), 
effect  of  species  and  temperature  gradients  in  the  flame  zone  on  the  spatial 
characteristics  of  the  plasma  (E/N,  electron  density  etc.) 

focus  on  plasma  radical  generation  in  the  preheat  zone  and  the  impact  on  overall  flame 
characteristics. 
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1 .  S.  Nagaraja,  V.  Yang,  I.  Adamovich,  “ Multi-Scale  Modeling  of  Pulsed  Nanosecond 
Dielectric  Barrier  Discharges  in  Plane-to-Plane  Geometry ,”  Journal  of  Physics  D: 
Applied  Physics  46  (15),  2013,  155205. 

2.  S.  Nagaraja,  V.  Yang,  Z.  Yin  and  I.  Adamovich,  “ Ignition  of  Hydrogen-Air  Mixtures 
using  Pulsed  Nanosecond  Dielectric  Barrier  Plasma  Discharges  in  Plane-to-Plane 
Geometry ,”  Combustion  and  Flame,  2013,  in  press 

3.  S.  Nagaraja,  and  V.  Yang,  “ Detailed  Comparison  between  Nanosecond  Plasma  and 
Thermal  Ignition  of  Hydrogen-Air  Mixtures”  to  be  submitted  to  Combustion  and  Flame. 

4.  S.  Nagaraja  and  V.  Yang,  “ Numerical  Investigation  of  Nanosecond  Plasma  Assisted 
Ignition  ofH2~,  CH4-  and  C2H4- Air  Mixtures”  to  be  submitted  to  Combustion  and  Flame. 

5.  S.  Nagaraja,  W.  Sun  and  V.  Yang,  “  Nanosecond  Plasma  Assisted  Ignition  of  n-Heptane- 
Air  Mixtures ,”  in  preparation. 
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Plasma-Assisted  Combustion 

Modeling  Framework 
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Model  Assumptions 


Governing  Equations 


•  Plasma  fluid  with  drift-diffusion 
approximation. 

•  Two  temperature  model:  electrons  at  Te 
(defined  using  mean  energy);  ions  and 
neutrals  at  gas  temperature,  Tgas 

•  Lookup  table  for  electron  transport  and 
rates  using  two-term  expansion  for 
electron  Boltzmann  equation 
(BOLSIG). 

•  Solution  to  mean-energy  equation  to 
update  electron  coefficients  at  each  time 
step. 

•  Uniform  pre-ionization  in  the  discharge 
volume.  No  photo-ionization  source 
term. 
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Validity  of  the  BOLSIG  approach  to  calculate  electron  rate  coefficients,  among  other  assumptions,  has 
been  validated  through  comparison  of  species  density  (O  and  OH),  temperature  and  input  energy  with 
experiments 
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Nanosecond  Plasma  Assisted  Ignition  and  Combustion 

Multi-Scale  Modeling  Framework 


School  of  Aerospace  Engineering 


PLASMA  DISCHARGE  DYNAMICS  FLOW  AND  COMBUSTION  DYNAMICS 


Plasma  Chemistry  Ignition  and  Combustion  Chemistry 

Ionization,  Excitation,  Dissociation,  Radical  initiation,  Chain  propagation, 


^  $  * 


•  ionization  wave 

•  quenching  of  excited 

•  cumulative  effects  of  multiple 

propagation 

species 

discharge  pulses 

•  electrical  breakdown 

•  ion  recombination 

•  convective  and  diffusive 

•  cathode  sheath  formation 

•  gas  heating 

transport 

•  electron  impact  dynamics 

•  ignition  and  combustion 

Solution  Algorithm 


Normalized  CPU  Time 
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Strategies  for  Computational  Efficiency 
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nanosecond 
voltage  pulses 


Adaptive  Time-Stepping 


•  At  varies  between  10'13  -10'12  s 

•  Semi-implicit  treatment  of  the  Poisson 
equation  to  circumvent  the  stiffness  arising 
from  tight  coupling  between  electric  field 
and  electron  density. 


•  At  fixed  at  10"9  s  . 

•  Electron  energy  equation  and  Poisson  equation  are 
not  solved  since  electric  field  effects  become 
negligible  and  the  space  charge  density  rapidly  decay 
as  the  applied  voltage  ends. 


Pulsing  Frequency  (kHz) 


Multi  Time-Scale  Treatment  of  Chemical  Source  Term 


•  speedup  by  40%  is  seen  by  using  the  multi  time-scale 
treatment  of  chemical  source  terms,  but  not  orders  of 
magnitude  speedup  observed  in  combustion  simulations 
without  plasma  discharge. 

•  at  high  pulsing  frequencies,  the  savings  with  using  HMTS 
reduce  because  more  time  spent  in  simulating  electric  field 
transients  during  breakdown  in  each  voltage  pulse. 
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*  Gou  et  al.,  Combustion  and  Flame  157  (2010)  1111-1121 
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Strategies  for  Computational  Efficiency 

plasma  combustion  chemistry  optimization 
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full  chem:  35  species,  287  reactions  optimized  chem:  19  species  111  reactions 


•  simulated  a  burst  of  8  nanosec  pulses  with  detailed  mechanism. 

•  removed  species  with  peak  mole-fraction  less  than  1  O'8 

•  ensured  that  E/N,  electron  and  radical  species  densities  and  temperature  (in  both  space  and 
time,  all  within  10%)  are  accurately  predicted  by  the  reduced  mechanism. 

•  provides  a  speed-up  of  ~  4  times  with  H2-air  plasma  ignition. 

•  expect  greater  savings  with  large  CxHy  mechanisms  and  2D/3D  simulations. 
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OSU  Plasma  Flow  Reactor 

40  -  160  torr,  300  -  500  K,  1  -  100  kHz 
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measurements 


ICCD  images  of  discharge  structure 


•  time  resolved  O  density  using  TALIF 
(uncertainty  +/-  30%) 

•  time  resolved  NO  density  using  LIF 
(uncertainty  +/-  30%) 

•  time  resolved  OH  density  using  LIF 
(uncertainty  +/-  20%) 

•  time  resolved  temperature  using 
rotational  CARS  and/or  LIF  thermometry 

•  ignition  delay  time  from  OH*  emission 
rise 

•  ICCD  imaging  of  discharge  structure  and 
flame  kernel  evolution. 
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Air  :  373  K,  60  torr,  40  kHz 
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OSU  Plasma-Coupled  Premixed  Flat  Flame 
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“direct  coupled”  configuration 


“plasma  upstream”  configuration 


•  low  pressure  ID  flame  (20  -  30  torr) 

•  H2/02/N2,  CH4/02/N2  and  C2H4/02/N2  premixed  flames 

•  FID  pulser:  14  kV  peak  voltage,  7  ns  FWHM,  ~3  mJ/pulse 

Measurements 

•  spatially  resolved  OH  density  using  LIF 

•  spatially  resolved  temperature  using  five-line  OH 
thermometry. 
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Nanosecond  Pulsed  Dielectric  Plasma  Simulations  in  N2/Air 
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Objectives 

•  Self-consistent  simulations  of  pulsed  nanosecond  discharges  in  air  with  detailed 
kinetics. 

•  Validation  with  experiments  and  analytical  model  results. 

•  Direct  insight  into  plasma  heating  and  radical  production  of  critical  importance 
in  combustion  applications. 


Model  Geometry 


Operating  Conditions: 


Applied  Waveform 


Left  Right 

Electrode  Electrode 


Initial  Pressure:  60  ton- 

initial  Temperature:  300  K 

Pulsing  Frequency:  40  kHz 

Gap  width:  1  cm 

Initial  Electron  Density:  107  cm'3 

Dielectric  thickness:  1.75  mm 

Dielectric  Constant:  4.3 

Pulse  Duration:  100  ns,  FWHM:  12  ns 

Peak  Voltage:  -22.5  kV  and  +17.5  kV 


S.  Nagaraja,  V.  Yang,  I.  Adamovich,  “Multi-Scale  Modeling  of  Pulsed  Nanosecond  Dielectric  Barrier 
Discharges  in  Plane-to -Plane  Geometry ,”  Journal  of  Physics  D:  Applied  Physics  46  (15),  2013,  155205. 
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Experimental  Validation  (Air  Discharge) 
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•  coupled  energy,  temperature  and  O  atom  density  predicted 
by  the  ID  model  are  within  20,  5  and  10%  of  experimental 
data,  respectively. 

•  coupled  energy  remains  fairly  constant  with  pulse  number, 
increasing  linearly  with  pressure,  and  nearly  independent 
of  pulsing  rates. 

•  O  atom  production  via  electron  impact  dissociation  and 
quenching  of  excited  N2  by  02  is  captured  accurately 
along  with  subsequent  decay  via  formation  of  03  over  ms 
timescales. 
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Detailed  Physics  over  ns-ms  Timescales 
air  discharge  (60  torr,  300  K,  40  kHz,  100  pulses) 
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temperature  evolution 
for  100  nanosecond  pulses 
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spatial  distribution  of  O  atom  density 
for  100  nanosecond  pulses 
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•  plasma  heating  effect  is  about  0.5  -  lK/pulse  in  air  and  nearly  independent  of  pulsing  frequency  (as  a  function  of 
pulse  number). 

•  rapid  gas  heating  produces  weak  acoustic  waves  which  propagate  into  the  gas  volume  from  both  ends.  The  strength  of 
these  waves  becomes  weak  as  overall  temperature  rises  from  heat  release  from  quenching  of  excited  species. 

•  a  fairly  uniform  temperature  profile  develops  in  the  plasma  volume  after  several  discharge  pulses,  owing  to  slow  but 
steady  (~0.5  K/pulse)  heat  release  primarily  from  relaxation  of  excited  species. 

•  repetitive  pulsing  results  in  efficient  production  of  atomic  oxygen  through  electron  impact  dissociation  during 
discharge  pulses,  and  quenching  of  excited  nitrogen  species  by  oxygen. 

•  volumetric  radical  generation  and  heating  by  pulsed  discharges  are  of  great  significance  for  ignition  and 
flame  stabilization  purposes. 
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Nanosecond  Pulsed  Dielectric  Plasma  in  H2-Air  Mixtures 
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operating  conditions 

Pressure:  40  -  160  torr 
Temperature:  373  -  600  K 
Pulsing  Frequency:  10-40  kHz 
Gap  width:  1  cm 
Initial  Electron  Density:  107  cm'3 
Dielectric  thickness  (Quartz):  1.75 
mm 

Dielectric  thickness  (Kalrez):  1.58  mm 
Dielectric  Constant  (Quartz):  4.3 
Dielectric  Constant  (Kalrez):  4-9 


applied  waveforms 


voltage  and  current 


Objectives 

•  comparison  of  OH  density  with  measurements  after  a 
burst  of  50  pulses. 

•  assess  the  accuracy  of  kinetics  model  at  low 
temperature,  pre-ignition  conditions 

•  detailed  investigations  of  NS  plasma  ignition  physics  and 
chemistry. 

•  sensitivity  of  ignition  process  to  key  system  parameters 
and  material  properties. 


S.  Nagaraja,  V.  Yang,  Z.  Yin  and  I.  Adamovich,  “ Ignition  of  Hydrogen-Air  Mixtures  using  Pulsed  Nanosecond  Dielectric 
Barrier  Plasma  Discharges  in  Plane-to-Plane  Geometry ,”  Combustion  and  Flame  (accepted). 
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H2-Air  Plasma  Combustion  Kinetics 
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nonequilibrium 
plasma  chemistry 

-► 

low  temperature 
radical  chemistry 

-► 

high  temperature 
combustion  chemistry 

•  low  temperature  (500-1000  K)  uncertainties  in  many  key  chain  branching  reactions. 

•  detailed  chemistry  mechanism  has  35  species  and  287  reactions. 

•  reduced  chemistry  mechanism  has  19  species  and  111  reactions. 
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E/N  and  Electron  Density  Evolution 
Pt  =  80  torr,  Tt  =  500  K,/=  60  kHz,  0  =  1.0,  FID  Pulser 
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•  breakdown  voltage  at  these  conditions  occurs  at  ~10  kV 

•  sharp  spike  in  current  is  seen  at  breakdown  before  it  drops 
rapidly  from  the  plasma  shielding. 

•  plasma  excited  species  production  happens  only  during  a  short 
duration  of  ~5  ns  when  E/N  is  high. 

•  E/N  and  electron  density  reach  a  periodic  steady  state  after  ~8 
pulses. 

•  nanosecond  discharge  efficiently  generates  radicals  and  excited 
species  during  each  pulse  because  of  high  peak  E/N. 
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key  low  temperature  pathways  for 
consumption  of  O  and  recirculation 
of  H  and  OH 


H  +  O.+M  -+HO,+M  O  +  HO.^OH  +  O, 


H  +  HO^OH  +  OH 


OH  +  O^H  +  O, 


model  predictions  for  OH  are  within  10%  of  measurements  in  H2-air  mixtures,  including  both  peak 
value  and  decay  rates. 

O  production  is  highly  sensitive  to  changes  in  eq.  ratio,  increasing  by  -50%  when  (j)  is  decreased 
from  0.12  to  0.06. 

H  and  OH  are  relatively  insensitive,  changing  by  -10  %  when  eq.  ratio  is  doubled. 
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How  is  ignition  achieved  with  nanosecond  plasma? 

Pt  =  80  torr,  Tt  =  500  K ,/=  60  kHz,  0  =  1.0,  FID  pulser,  115  pulses  ~  2  ms 
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temperature  evolution  vs  time 
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spatial  evolution  of 
OH  and  temperature 
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•  heat  transport  plays  a  minor  role,  local  plasma  chemistry  effects  are  critical 
in  producing  “volumetric”  ignition 

•  secondary  peaks  in  OH  density  near  the  boundaries  is  generated  from  H02 
which  has  accumulated  due  to  low  temperatures 


4^  Ol  Ov  3 

OH  Number  Density  (cm' ) 


Number  Density  (10  cm 3)  Number  Density  (10  15  cm' 
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Spatial  Evolution  of  Radicals  during  H2-Air  Ignition 

Pi  =  80  torr,  Ti  =  500  K ,/=  60  kHz,  0  =  1.0,  FID  pulser,  115  pulses  ~  2  ms 
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O  H  OH 


ho2 


•  a  small  increase  in  temperature  near  ignition  significantly  increases 
the  chain  branching  reaction  rates. 

•  radical  concentration  profiles  are  much  steeper  than  the  temperature 
distribution,  with  well  pronounced  maxima  near  the  centerline. 

•  both  O  and  H  densities  increase  by  ~3  times  within  0.4  ms  near 
ignition,  with  OH  density  increasing  by  4  times. 

•  low  temperatures  at  the  boundaries  because  of  heat  losses  aid  the 
accumulation  of  H02  which  generates  OH.  The  secondary  peaks  in 
OH  profiles  near  the  boundaries  result  from  this  pathway. 
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H2-air  pulsed  nanosecond  plasma  ignition 

what  can  we  infer  from  emission  images? 
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OSU  Experiment 

Pi  =  104  torr,  T}  =  473  K ,/=  40  kHz,  <t>  =  1 .0,  CPT  pulser 


l;t:smc.  5  ms  l-'himg,  5,1  inn  flume,  5,4  ms 


ignition  first  observed  near  edges,  within  0.1  ms,  ignition  observed  the  entire  volume  ignited 
due  to  higher  electric  fields  all  along  the  centerline.  within  0.4  ms  of  first  detected 

flame  emission. 


•  the  present  ID  model  simulates  a  particular  cross-section. 

•  although  the  ID  model  cannot  capture  edge  effects,  it  is  able  to  explain  the  spreading  of  the  ignition 
kernel  from  the  centerline  towards  the  boundaries. 

•  predictions  are  in  line  with  observations  that  local  plasma  chemistry  dominate  over  heat  transport  effects 
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Nanosecond  Plasma  Ignition  vs  Thermal  Ignition 

is  there  any  difference? 
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species  evolution  at  center  for 
NS  ignition 


ignition  delay  times  with  nanosecond 
plasma  and  a  volume  heat  source 


Pt  =  80  torr 
7)  =  500  K 
/=  60  kHz 
0=  1.0 
FID  pulser 
115  pulses 
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steps  in  the  nanosecond  plasma  ignition  process 


N2  +  e  N2(A3,  B3,C3)  N2(A3)  +  02  ^  O  +  O(D)  O  +  OH  H  +  02 

H2  +  e^H  +  H  ^  N2(B3,C3)  +  H2  ^  H  +  H  ^  0  +  H02^0H  +  02  ^  H2  +  0H^H20  +  H 

02  +  e^0  +  0  H  +  H02  ^  OH  +  OH 


•  the  high  activation  energy  chain  initiation  reactions  are  replaced  by  electron  impact  reactions 
with  NS  plasma. 

•  for  the  same  input  energy,  thermal  ignition  delay  is  ~  60%  higher. 

•  plasma  generated  radicals  trigger  heat  release  from  fuel  oxidation  at  -700  K,  as  opposed  to  auto¬ 
ignition  temperature  of  -960  K  under  same  conditions. 
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Effect  of  Burst  Size  and  Dielectric  constant  on  Ignition 


iv. 
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ignition  delay  vs  #  pulses  in  burst 

Pi  =  80  torr,  7J  =  473  K 
/=  40  kHz,  0=1.0,  CPT  pulser 


ignition  delay  sensitivity  to  eq.  ratio 

Pt  =  80  torr,  Tt  =  500  K 
/=  60  kHz,  FID  pulser 


•  there  is  a  minimum  #  of  pulses  in  burst,  below  which  no  ignition  is  observed. 

•  ignition  characteristics  are  highly  sensitive  to  dielectric  properties. 

•  uncertainty  in  the  dielectric  constant  values  should  be  considered  during  the  validation  process. 

•  ignition  delay  reduction  with  increase  in  burst  size  is  especially  pronounced  for  lean  mixtures. 
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Effect  of  Pressure  and  Pulsing  Rates  on  H2-Air  Ignition 
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Tt  =  473  K ,/=  40  kHz,  <|)  =  1.0,  CPT  pulser 
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•  ignition  delay  reduction  with  increase  in  pressure  is  well  reproduced  by  the  model. 

•  increase  in  pressure  results  in  nearly  linear  rise  in  input  energy  per  pulse  because  of  its  dependence  on  number  density. 
Faster  addition  of  energy  results  in  more  rapid  ignition  at  higher  pressures. 

•  the  nonlinear  trend  of  #  pulses  required  for  ignition  as  a  function  of  pulsing  frequency  is  not  reproduced  by  the  model. 

•  the  model  predicts  that  input  energy  per  pulse  is  nearly  independent  of  pulsing  frequency,  which  may  not  be  tme. 

•  lowering  of  input  energy,  because  of  residual  electron  density  effects,  with  rise  in  pulsing  rates  may  explain  the 
observed  nonlinear  trend. 
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OH  Density  Decay  after  50  Pulse  Burst 
in  CH4-,  and  C2H4-Air  Mixtures 
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Pt  =  100  torr,  Tt  =  500  K ,/=  10  kHz,  FID  pulser 

CH4-air 


C2H4-air 


CH4-air 
(64  species) 

GRI  Mech  3.0 
CH4/N2/02  plasma 
NOX  reactions 

C2H4-air 
(70  species) 

use  Mech 
C2H4/N2/02  plasma 
NOX  reactions 


•  GRI  Mech  3.0  has  been  validated  extensively  in  1000-2500  K  and  25  torr  to  10  atm  range. 

•  USC  Mech  has  been  validated  in  900-2500  K  and  16  torr  to  10  atm  range. 

•  model  consistently  over-predicts  OH  density  by  -50%  in  CH4-air  mixtures. 

•  growth  rate  is  correctly  predicted  in  C2H4-air  mixtures,  but  the  decay  rate  in  slower  than  measurements. 

•  low  temperature  uncertainty  in  chain  reactions  may  be  the  primary  reason  for  deviations 


Ongoing  Work 

•  different  CH4-  and  C2H4-air  combustion  chemistry  integrated  with  plasma  kinetics  are  being  tested  to 
assess  their  relative  performance  on  predicting  low  temperature  radical  production/decay. 


Temperature 
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Nanosecond  Plasma  Ignition  of  nHeptane-Air 

(in  collaboration  with  Wenting  Sun) 


School  of  Aerospace  Engineering 


Objective 


i 

r 

•  understand  the  effect  of  NS  plasma  on  n- 

1st  stage 

f  2nd  stage 

heptane-air  ignition  chemistry  through 

R-H  R  (H  abstraction) 

1  H202  -»  20H 

self-consistent  simulations. 

R  +  02  R02  (exothermic) 

/  CH20  -»  CO,  H20 

1/  CO  ->  co2 

•  investigate  the  effect  of  radical  addition  to 

R02  -»  H02,  H202,  CH20  etc. 

) - r 

the  “low  temperature”  and  “high 
temperature”  steps  of  the  2-stage  ignition 
process. 

time 


nC7H16-air  plasma  combustion  kinetics  (154  species) 


C7H16/N2/02  combustion 
(LLNL  reduced  mech  +  NOX  reactions) 
n  -  C7H]6  +OH^>  2-  C7Hl5  +  H20 
C7Hl4OOH  <->  C7HuO  +  OH 

c2h3  +o2  ->  ch2cho + o 

CH20  +  OH  <->  HCO  +  H20 

h+o2^o+oh 

H202+M  -+OH  +  OH  +  M 


C7H16/N2/02  plasma  reactions 


e  +  N2^N2(A\B\C\al)  +  e 
e  +  N2^N  +  N(2D)  +  e 
e  +  02  — ^  O  +  O(-D)  +  e 
e  +  C-jH 16  — ^  Cy//j5  +  H  +  e 
e  +  C7Hl6  — ^  C^H33  +  CH3  +  e 
e  +  C-,H,  r  — ^  CrH, ,  +  C-.H .  +  6 


N2{A3)  +  02  -^n2+o  +  o 
N2  (A3 )  +  C7Hl6  -^N2  +  C7Hl5  +  H 
N2  (A3  )  +  C7H16^N2+  C7Hl4  +  H2 
N2  (A3 )  +  C7Hl6  -^N2  +  C6Hl3  +  CH3 
0(D)  +  C7Hi6  -»  C7Hl5  +  OH 


*  C7H16  (electron  impact  and  with  excited  species)  reaction 
rates  estimated  from  C3H8  based  plasma  reactions. 
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Nanosecond  Plasma  Ignition  of  nHeptane-Air 

Effect  of  NS  Pulses  on  1st  Stage  Delay  Time 
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model  geometry 

discharge 

gap 


operating  conditions 


O,  H  and  OH  density  evolution 
at  center  of  discharge  volume 
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temperature  evolution 
at  center  of  discharge  volume 


“self-acceleration”  of  low  temperature  chemistry 


addition  of  small  amount 
of  radicals  accelerates  the 
H  abstraction  step 
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R02  produces  more 
radicals  which  accelerate 
the  whole  process  further. 
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Nanosecond  Plasma  Ignition  of  nHeptane-Air 

Effect  of  NS  Pulses  on  Overall  Ignition  Delay  Time 


“staggered”  application  of  NS  pulses 


mn 


NS  pulser  off 
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operating  conditions 

Pi  =160  torr 
Ti  =  600  K 
/=  60  kHz 
0=  1.0 

8  kV  Gaussian  pulses 
10  ns  duration 


only  a  few  NS  pulses 
sufficient  to  rapidly  trigger 
1  st  stage  temperature  rise 


25  NS  pulses  are  applied  after 
the  1  st  stage  to  reduce  the 
overall  ignition  delay 
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•  the  “staggered”  application  of  NS  pulses  result  in  ~  40% 
reduction  in  ignition  delay  time 

•  it  is  evident  that  the  2nd  stage  is  less  sensitive  to  radical 
addition  by  NS  pulses  than  the  1st  stage. 

•  heating  provided  by  the  NS  pulses  after  the  1  st  stage 
accelerate  the  decomposition  of  H202  and  reduce  ignition 
delay. 

H202  +  M->  OH  +  M 

•  inclusion  of  NOX  catalytic  reactions  change  the  predictions 
by  ~5%  because  of  following  new  OH  generation  pathways. 

NO  +  H02->  N02  +  OH 

NO  +  CH302->  N02  +  CH20  +  OH 
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Time  (s) 
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Mole  Fraction 
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Nanosecond  plasma  coupled  premixed  flame 
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CH4-air 
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physical  setup 

Wiremesh 

1 1 1 1 1 1 1 1 1 
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HVE 

Flame 

Pulse 

Discharge 

Generator 


t  t  t  t  t  t  ft  t  t 

Premixed 


Mixture 


operating  conditions 

Pressure:  25  torr 

Inlet  Temperature:  650  K 

Eq.  ratio:  1.07 

Gap  width:  4.0  cm 

initial  Electron  Density:  1 07  cm'3 

-15  kV  peak  voltage 

7  ns  FWHM 

Mdot :  0.00377  kg/m2-s 


validation  of  flame  model 
with  CHEMKIN  solution 


Height  above  Burner  (cm) 


CH4-air  plasma  flame  kinetics  (75  species) 

GRI  Mech  3.0  +  CH4/N2/O2/CO/CO2  plasma  +  NOX  chemistry 

•  electron  impact  processes  of  both  reactant  (CH4,  02,  N2) 
and  product  species  (H2Q,  CO,  C02)  considered. 


N2+e^N2(A2,B\C3,al)  +  e 
02  +  e  — ^  0  +  0  +  6 
CO  +  e—>  CO*+e 
H2  0  +  6  — ^  O  +  H2 


what  are  most  important 
plasma  pathways  pertaining  to 
H20,  CO,  C02  ? 
is  plasma  species  production 
in  preheat  zone  more 
important  than  downstream? 
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Plasma  Coupled  Premixed  CH4-Air  Flame 

Pt  =  25  torr,  Tt  =  650  K,  @=  1.07 
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physical  setup  reduced  electric  field  vs  height 
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electron  density  vs  height 


27.5 
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22.5 
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0.5  1  1.5 

Height  above  Burner  (cm) 


high  E/N  (100  -  600  Td)  and  high  electron 
densities  (~lel3  cm-3)  allow  for  efficient  radical 
generation  by  NS  pulses  in  preheat  zone,  where 
they  may  provide  significant  benefit. 


note  that  radicals  (O,  H,  OH  etc)  concentration  in 
the  flame  zone  is  already  of  the  order  10,000  ppm, 
so  plasma  cannot  make  much  impact 


•  high  E/N  downstream  of  the  flame  can  be  attributed  to  high  temperatures  and  low  number  density. 

•  E/N  in  the  preheat  zone  reaches  ~  600  Td  at  22.5  ns.  Plasma  radical  generation  in  this  zone  may  have  a  significant 
impact  on  flame  characteristics. 

•  sharp  peak  in  the  E/N  profile  at  20  ns  at  the  right  boundary  indicates  the  cathode  sheath  region. 

•  electron  density  distribution  is  fairly  uniform  in  the  entire  domain  reaching  peak  value  of  2xl013  cm-3  at  27.5  ns. 

•  total  input  energy  during  the  pulse  was  2.7  mJ 


2  cm 
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Plasma  Coupled  Premixed  CH4  -  Air  Flame 

P:  =  25  torr,  Tt  =  650  K,  @=  1.07 
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physical  setup 
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N2  excited  species  vs  height 


O  atom  density  vs  height 


plasma  generation 
of  O  atoms  in  the 
preheat  zone  may 
provide  significant 
benefits. 


Height  above  Burner  (cm) 
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•  the  production  rates  of  N2(A3)  and  N2(B3)  in  the  preheat  zone  is  about  2  times  higher  than  downstream 
because  of  higher  N2  number  density. 

•  electron  impact  dissociation  of  02  in  the  preheat  zone  results  in  ~30  times  increase  in  O  atom  density  within 
30  ns 

•  the  excited  species  are  quenched  rapidly  after  the  pulse  resulting  in  further  production  of  O  and  other  radicals 

ongoing  work 

•  we  are  performing  longer  timescale  simulations  to  understand  the  effect  of  repetitive  application  of  discharge 
pulses  on  flame  dynamics. 

•  the  effect  of  NS  discharges  on  H2-air,  CH4-air  and  C2H4-air  premixed  flames  are  being  investigated. 

•  close  collaboration  with  OSU  group  is  pursued  for  obtaining  greater  insight  through  experiments  and  high  fidelity 
modeling 
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Where  we  go  from  here? 
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High  fidelity  ID  numerical  tools  for  construction  and  validation  of 
robust  plasma  combustion  kinetic  models 

•  detailed  studies  of  the  plasma  coupled  premixed  flame  system  for  a  variety  of  fuels. 

•  development  of  the  counterflow  plasma  flame  simulation  framework. 

•  close  collaboration  with  other  MURI  team  members  for  model  validation  and  critical 
assessment  of  the  plasma  combustion  kinetic  models 

2D/3D  simulations  of  nonequilibrium  plasma  in  complex  flow  environments 

•  High  fidelity  simulations  of  single  filament  discharge  in  2D  with  detailed  chemistry. 

•  Large  Eddy  Simulation  (LES)  of  H2  jet  in  supersonic  02  crossflow  in  the  presence  of  a 
nanosecond  plasma  source. 

•  theoretical  framework  to  understand  plasma-flow  interactions. 


